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Abstract. The magnetization distribution, its energetic characterization by the interlayer coupling con- 
stants and lattice dynamics of (OOl)-oriented Fe/Pt multilayers are investigated using density functional 
theory combined with the direct method to determine phonon frequencies. It is found that ferromagnetic 
order between consecutive Fe layers is favoured, with the enhanced magnetic moments at the interface. The 
bilinear and biquadratic coupling coefficients between Fe layers are shown to saturate fast with increasing 
thickness of nonmagnetic Pt layers which separate them. The phonon calculations demonstrate a rather 
strong dependence of partial iron phonon densities of states on the actual position of Fe monolayer in the 
multilayer structure. 
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The technological importance of metallic multilayers con- 
sisting of alternating layers of magnetic and nonmagnetic 
material stimulated both experimental and theoretical stud- 
ies of these systems which continues over the last two 
decades. The properties of magnetic multilayers are to 
a large extent determined by the interface properties. In 
order to understand whether future applications of mul- 
tilayers as magnetic data storage materials are possible, 
their microscopic properties have to be investigated to de- 
termine the magnetic interactions on the nanoscale and 
to understand the dependence of these interactions on the 
multilayer structure. Since the discovery of composition 
and annealing dependence of the magnetic properties in 
Fe/Pt bilayers [1 , a number of experimental papers have 
been devoted to the magnetic properties of the Fe/Pt mul- 
tilayers [2113]. 

The materials themselves are of great technological 
importance since they are characterized by large uniax- 
ial magnetocrystaline anisotropy energy |4j , which in turn 
leads to high temperature stability of the data magnet- 
ically recorded in such a medium [5]. The magnetocrys- 
taline anisotropy energy density has been already anal- 
ysed [313] using ab initio techniques. However, to the best 
of our knowledge, the dependence of the type of magnetic 
ordering on the thickness of Fe and Pt layers, as well as 
lattice dynamics in these systems, have not been investi- 
gated with modern computational techniques until now. 



e-mail: Pawel.Jochym@ifj.edu.pl 



The purpose of the present paper is to analyse mag- 
netic ordering and lattice dynamics in the (OOl)-oriented 
Fe/Pt multilayers using an ab initio approach. We have 
computed the magnetization distribution in a series of rep- 
resentative multilayer structures with one up to twelve 
atomic monolayers (MLs) in an iron layer, and one up to 
five platinum MLs in the spacer layer. The calculations 
were performed self-consistently and the configurations 
with ferromagnetic (FM) and antiferromagnetic (AF) or- 
der between the consecutive iron layers were obtained and 
their energies and magnetization distributions were com- 
pared with each other. We also report the values of bilinear 
and biquadratic coupling constants, which quantify the 
energy difference between the AF and FM configuration 
of the Fe layers within the model introduced by Fuller- 
ton and Bader [7^, and depend on the thickness of both 
magnetic iron layers and (almost) nonmagnetic platinum 
layers. 



The paper is organized as follows. The computational 
methods are shortly described in Sec. II. In Sec. Ill we 
present the magnetic properties of Fe/Pt multilayers, and 
systematically investigate the energy difference between 
the AF and FM configurations, depending on the num- 
ber of Fe and Pt atomic MLs in the considered multilayer 
structure. The phonon dispersions and the phonon densi- 
ties of states for two representative Fe/Pt multilayer sys- 
tems are analysed in Sec. IV. We present a short summary 
and the main conclusions in Sec. V. 
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Fig. 1. Side view of a representative 5/5 Fe/Pt multilayer 
structure before geometry optimization (with five atomic MLs 
of iron and five atomic MLs of platinum) . Non-equivalent MLs 
in the unit cell are labeled by characters A — F; the unit cell 
in the FM structure, used also for phonon calculations, is indi- 
cated by dashed lines (the magnetic unit cell in an AF config- 
uration is twice larger) . The bottom row of Pt atoms is shown 
only for clarity. 



2 Computational methods 

The calculations were carried out with standard density 
functional theory (DFT) code as implemented in the VASP 
package [S]. We have used well established ultrasoft gen- 
eralized gradient approximation (GGA) pseudopotentials 
[5] with Perdew-Wang 91 gradient correction scheme (as 
given in Ref. [TU]) as input for phonon calculations. The 
calculations of magnetic structure of the system were car- 
ried out with pseudopotentials introduced by Perdcw, Burke 
and Ernzerhof [11,12] (so-called PBE-GGA pseudopoten- 
tials) as implemented in the VASP package [9]. The differ- 
ence in method used for these two calculations does not 
influence our general conclusions since the results from dif- 
ferent methods are not compared with each other. The cal- 
culations neglected spin-orbit coupling which may some- 
what influence the results, especially for Pt atoms, since 
they possess fairly strong spin-orbit component. 

The series of structures investigated in this paper was 
derived from the basic (001) bcc iron lattice with plat- 
inum atoms replacing iron atoms within platinum layer 
(by layer we mean one or more MLs of the same material) . 
In this way some MLs of Fe atoms were replaced by Pt 
ones. Then the entire multilayer structure was optimized 
with respect to atomic positions and its unit cell size. The 
convergence criterion has been set at root mean square 
(RMS) force < lO"'^ eV/A and RMS stress < 10 MPa. The 
atomic forces and magnetic moment density distribution 
have been calculated using final optimized structures. One 
example of such a structure before optimization is shown 



in Fig. [TJ its unit cell in the FM (ground state) config- 
uration contains atoms belonging to five atomic MLs of 
iron and to five atomic MLs of platinum, in so-called 5/5 
system. 

The structure optimization led to substantial changes 
in the structure of the considered multilayer systems. Namely, 
the inter-layer distance in the Pt layer increased to the 
point where this layer adopts the well known experimen- 
tally fee structure. This change is illustrated in Fig. ^ 
showing the local c/a ratio of lattice constants for each 
atomic monolayer in the structure. It is straightforward 
to see by analysing these data that the structure of the 
platinum layer changes from the starting bcc structure 
(assumed at the beginning of the optimization process) 
to the fee lattice - well-known experimentally. This may 
be interpreted, together with phonon calculations below, 
as a strong evidence of rather realistic description of the 
physical properties of the multilayer systems in the pre- 
sented calculations. The above change starts to show up 
for quite thin Pt layers — even the system consisting of 
four Pt atomic MLs shows already the interlayer spacing 
close to the fee structure. For thicker Pt layers this effect 
is very strong (see Fig. [2|) . The above finding agrees quite 
well with experimental findings of Antel et al. [13] , where 
a similar process of gradual change from tetragonally dis- 
torted bcc (bet) to tetragonally distorted fee (fet) lattice 
was observed. The iron layer is quite stable but when it is 
thin (few atomic MLs) it is prone to be distorted towards 
the fee structure. This tendency may be also observed in 
Fig. [21 where the distortion of c/a ratio goes deep into 
the Fe layer, while almost the whole Pt layer keeps its 
optimum c/a ratio. Similar structural changes in Fe/Pt 
systems have also been observed experimentally [T31I1]. 

All calculations were carried out for spin-polarized sys- 
tems with the self-consistent optimization of the magne- 
tization distribution within a given Fe/Pt multilayer sys- 
tem. The set of calculated structures contains a range of 
variants from the simplest 1/1 structure (periodic struc- 
ture with 1 ML of Fe and 1 ML of Pt in its unit cell) up 
to 4/8 (i.e., with 4 atomic MLs of Fe and 8 atomic MLs 
of Pt) in increments of 2 MLs, and also, in addition, 1/12 
structure (1 ML of Fe and 12 MLs of Pt). 

The magnetic properties of Fe/Pt multilayers were in- 
vestigated systematically (see Sec. Ill) by changing the 
thickness of Fe and Pt layers in a periodic system con- 
sisting of four layers, i.e., twice larger than that shown in 
Fig. [l] with two Fe and two Pt spacer layers interchanged 
with each other, with periodic boundary conditions. Sev- 
eral different starting configurations for the magnetization 
distribution were used in all cases in order to obtain un- 
biased results for a given system. Here we have used a 
similar approach to that employed to investigate Fe/FeSi 
multilayers in our recent paper |14| . 

Starting from equilibrium multilayer configurations ob- 
tained as a result of electronic structure calculations, we 
investigated lattice dynamics for two representative Fe/Pt 
multilayers: 4/4 and 5/5 systems (see Fig.[T]), reported in 
Sec. IV. Thereby we used the so-called direct method of 
lattice dynamics calculation [TS]. It follows the early idea 
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Fig. 2. Inter-layer distance in the 4/8 structure after struc- 
ture optimization. Dashed lines indicate c/a ratio for ideal bcc 
(bottom line) and fee (top line) structures. The shaded areas 
indicate the extension of the platinum layers. 



of calculating the force constants with ab initio methods 
|16| by implementing the full crystal symmetry given by 
the crystallographic space group [TS] . As such, it is able to 
achieve high accuracy and became a very powerful compu- 
tational tool in the last decade. By construction it allows 
one to determine accurate phonon frequencies at the high 
symmetry points in the reciprocal space, and gives other- 
wise the phonon spectra with high accuracy, if the force 
constants decrease fast enough with increasing distance, 
as in the present case. 

The formulation of Ref. [15] provides an algorithm to 
calculate force constants automatically for a system of any 
size. The procedure employed in the present study used 
displacements of 0.03 A in each non-equivalent direction 
for every non-equivalent atom. The number of necessary 
displacements was thus determined by the system symme- 
try and by the size of the multilayer unit cell. Next, the 
matrices of force constants were found using the singu- 
lar value decomposition method. Finally, the dynamical 
matrix was constructed and diagonalized for each wave 
vector k. 



3 Magnetic properties 

When Fe layers are separated by layers of nonmagnetic Pt 
spacer in a Fe/Pt multilayer, the first question concerns 
the distribution of magnetization over the Fe atoms. In 
the case of platinum spacer we have found that the inter- 
face proximity substantially enhances magnetic moments 
of iron atoms. This can be clearly observed for the 4/4 
Fe/Pt structure, with magnetic moments m ~ 2.8 fis and 
~ 2.5 hb for the Fe atoms at the interface and in the in- 
ner Fe ML (Fig. ^ , respectively, and with differences less 
than 2% between the FM and AF configuration. 

Also in larger systems we have found that the magnetic 
moments m at the interface are systematically enhanced 
and are rather close to 3.0 ^s, i.e., the 3d electrons are 
almost fully polarized. In contrast, the moments at the Fe 
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Fig. 3. Distribution of magnetic moments rrii in the 4/4 mul- 
tilayer structure consisting of two layers with four Fe MLs sep- 
arated by Pt spacer with four Pt MLs, as obtained for two 
magnetic configurations: (a) two magnetic unit cells in the FM 
ground state, and (b) magnetic unit cell of the (locally stable) 
AF order between two neighbouring Fe layers. 



atoms of inner MLs are smaller and gradually approach 
the bulk value, ttt-fc — 2.2 ^b, with increasing thickness 
of the Fe layer. While the later behavior may be expected 
and in fact is observed experimentally [T^], the magne- 
tization increase at the interface Fe atoms characterized 
by a lower coordination of Fe atoms is not obvious. We 
interpret it as a manifestation of electron correlation ef- 
fects which increase local moments at Fe atoms when the 
kinetic energy of 3d electrons is decreased by a lower Fe 
coordination |17) . In agreement with this interpretation, 
the moment enhancement at the interface becomes even 
more pronounced in the extreme case of a single iron ML, 
where Fe atoms have only platinum coordination and their 
magnetic moments are increased up to the level of satu- 
ration m ~ 3.0 /is, i.e., by 35% over the bulk iron value 
TTipo — 2.2 fiB (the magnetic moments of Pt atoms at the 
interface are here quite small mpt — 0.3 hb)- 

The above magnetic moment amplification effect at the 
interface in case of Pt spacer is remarkably different from 
the behaviour found recently for Fe/FcSi multilayers |14) . 
There, the magnetic polarization at interface Fe atoms is 
reduced, while the magnetic moments at Fe atoms of the 
second ML are somewhat enhanced. This demonstrates, 
that the moment enhancement may occur only in case 
of a metallic spacer, while we suggest that the electronic 
states at an insulating interface behave more like molec- 
ular bonding and thus 3d electrons are partly bound and 
the magnetization, is reduced. 

In addition, large magnetic moments at Fe atoms of 
the metallic interface induce weak FM polarization in the 
neighbouring Pt ML {m ~ 0.4 fj,B m case of 4/4 struc- 
ture of Fig. [3]) . The induced magnetization of the Pt layer 
has been observed experimentally \13\ to be on the same 
level as calculated here. There is however a difference in 
the distribution of the magnetic moment inside Pt layer 
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— the present results indicate large enhancement at the 
interface, while the experimental results suggest more uni- 
form distribution of the moment across the whole Pt layer. 
This difference may be connected with the contribution 
of spin-orbit interaction which, according to the experi- 
mental results [13], is absent or less pronounced in the 
case of iron-transition metal multilayers. It might be also 
to certain extent attributed to interdiffusion in the mea- 
sured sample, following the suggestion of Ref. [13] . Unfor- 
tunately, such effects cannot be included in the electronic 
structure calculations at present. 

Our main goal in investigating the magnetic proper- 
ties of multilayers is to establish the dependence of the 
favoured magnetic ordering in the underlying system on 
the thickness of the spacer (here Pt) layer and to char- 
acterize it by certain global parameters. Here we follow a 
phenomenological model of FuUerton and Bader [7] intro- 
duced to describe the energy difference between an AF and 
a FM configuration of a given multilayer structure. In this 
model it is assumed that the energy difference between an 
AF and a FM configuration of the system, 



AE = Eaf - Ei 



FM, 



(1) 



may be expressed by the following simple formula: 



E 



AF(FM) 



Ji mi • m2 + J2 (mi • mj)^ + Eo + E,„. (2) 



Here Ji and J2 are the bilinear and biquadratic coupling 
parameters, respectively, while mi and 1112 are magnetic 
moments of two interacting layers normalized per area 
unit (only the interactions between the nearest neighbour 
magnetic layers are considered). Eq and Em are reference 
energies of nonmagnetic structure and magnetic energy of 
the isolated layers, respectively. Note that the magnetic 
moments depend on the number of Fe MLs, and are prac- 
tically the same for the AF and the FM configuration of 
a given multilayer structure (see below). 

The derivation of the parameters {^1,0/2} in Eq- © 
requires the calculation of energies of the (locally stable) 
AF, FM and nonmagnetic variants of each investigated 
structure separately. Here we use the same method as the 
one introduced recently for Fe/FeSi multilayers ^14, . Note 
that the same result could also be obtained by calculations 
performed for two systems with spins oriented in two dif- 
ferent directions (e.g. the reference system and the one ro- 
tated by 90 degrees), but this would require employing the 
non-coUinear spins technique which is more involved and 
more technically demanding than the approach used here. 
In each case we determined these locally stable configu- 
rations and next analysed the results in order to extract 
information about the interlayer exchange couplings. 

Furthermore, the procedure used here requires the val- 
ues of magnetic moments m; of two coupled layers I — 1,2. 
The magnetic moments have been derived from the cal- 
culated spin-polarized electronic structure by numerical 
integration of the magnetic moment density in the spher- 
ical volume with radius r = 1.3 A around each atom. 
These individual moments were subsequently summed up 
within respective layers. We have verified that the result- 
ing moments for each layer were practically the same (the 
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Fig. 4. Energy difference AE ([T]) between AF and FM order 
of Fe layers in Fe/Pt multilayers for 1—4 atomic Fe MLs, as a 
function of increasing thickness of the Pt spacer layer. 



difference is less than < 1%) for both AF and FM con- 
figuration of the same structure. Further details concern- 
ing this procedure of calculating the coupling parameters 
{JitJ2} were described in the previous paper [14) . 

The energy difference AE ^ is presented in Fig. H] as 
a function of the increasing thickness of platinum spacer 
layer for several Fe/Pt multilayers. The data shows a clear 
preference towards the FM ordering of consecutive Fe lay- 
ers already for the systems with only up to four Fe MLs. 
Note, however, that the energy difference AE decreases 
(with the exception of a single Fe ML) with increasing 
thickness of the iron layer, while the platinum layer con- 
sists of a single ML. 

The above behaviour is again different from the Fe/FeSi 
system, where the order between the Fe layers is AF for 
thin layers, and becomes FM only for sufficiently thick 
FeSi spacer. Furthermore, it is quite remarkable that the 
energy difference AE for all the considered Fe/Pt mul- 
tilayers except for the system with one Fe ML converges 
quite fast to a constant value with the increasing thickness 
of the Pt layer. This feature appears to be again universal 
with exception of the system with a single iron ML which 
does not contain bulk-like iron atoms with only iron co- 
ordination, and thus is expected to behave in a different 
way. 

The thickness dependence of the coupling coefficients 
between Fe layers Ji and J2, see Eq. ^, are presented 
in Fig. [5] It is interesting to compare them with the ones 
obtained before for metal-insulator Fe/FeSi multilayers. 
There are some similarities and some differences between 
the characters of their thickness dependence in Fe/FeSi 
multilayers [13], and in the present metallic Fe/Pt ones. 
First, in both cases, one finds rather weak dependence of 
J2 on the thickness of the spacer layer (top right panel 
of Fig. [SJ . In case of one or two Fe MLs there is also no 
noticeable convergence of J2 with increasing thickness of 
the spacer layer, while for thicker Fe layers the values of 
J2 are very small indeed for 8 Pt MLs and suggest that 
this parameter becomes gradually irrelevant for increas- 
ing thickness of spacer layers. Second, there is a qualita- 
tive difference in J2 behaviour as a function of iron layer 
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Fig. 5. Bilinear Ji (bottom) and biquadratic J2 (top) coupling 
coefficients defined by Eq. as obtained for Fe/Pt multilayers 
for increasing thickness of iron layer (left) and platinum spacer 
layer (right). 



thickness with the one found in Fe/FeSi systems. In the 
case of Fe/FeSi multilayers one finds consistent and close 
to the exponential decay of J2, while for Fe/Pt we could 
spot some kind of oscillatory shape in the J2 plot (top left 
panel of Fig. [S]). This we interpret as a manifestation of 
the metallic behaviour, in contrast to the gradual decay 
of the coupling parameter J2 in an insulating medium. 

On the contrary, the behaviour of the bilinear coupling 
parameter Ji (bottom part of Fig. [5]) shows more similar- 
ities between Fe/FeSi and Fe/Pt case. In both systems 
the values of the Ji converge quite fast with increasing 
thickness of both iron and spacer layer. In case of the 
iron-thickness dependence (bottom left panel of Fig. [5]) 
the asymptotic value appears to be close to zero, while 
the data suggest that it is finite for the spacer-thickness 
dependence (bottom right panel of Fig. [SJ. 

A single ML of iron is again remarkably different from 
other multilayer systems with thicker Fe layers. The val- 
ues of J2 coupling constant are in this case negative and 
rather large, in contrast to thicker Fe layer. Furthermore, 
the values of | J2I do not decrease with increasing thick- 
ness of Pt spacer, unlike for thicker layers. Also the val- 
ues of Ji show a different behaviour from that found for 
ticker layers, and the change to positive Ji found at 8 Pt 
MLs followed by a large negative value for 12 Pt MLs is 
even somewhat surprising. It suggests that the magneti- 
zation density within Pt layers exhibits oscillations, which 
in some cases may even lead to weak AF coupling (see also 
Fig.U). Altogether, all these features confirm that the be- 
haviour of a single Fe ML is dominated by the interface, 
and follows from a rather strong FM polarization of Pt 
atoms at the metallic interface, which extends into the Pt 
layers. 




Fig. 6. Phonon dispersion relations of the Fe/Pt multilayer 
structures in FM configuration as functions of wave vector k 
along high symmetry directions in the Brillouin zone for: (a) 
4/4 multilayer, and (b) 5/5 multilayer of Fig.[T] The high sym- 
metry points in units of ^ are: F = (0,0,0), X = (1,0,0), 
M = (1,1,0),Z = (0,0,1)." 



4 Lattice dynamics 

We complete the discussion of the physical properties of 
Fe/Pt multilayers by presenting the phonon spectra for 
selected optimized structures. In order to reach unbiased 
conclusions, the lattice dynamics has been calculated for 
two systems, using the 2 a/2 x 2\/2 x 1 superccU for the 
4/4 multilayer structure, and the \/2x V2 x 1 supcrcell for 
the 5/5 one, respectively (the side view of the 5/5 struc- 
ture is shown in Fig. [Ij. These calculations are rather 
demanding and do not allow us to treat more atoms in 
a unit cell at present, but the unit cells of these two sys- 
tems are already large enough to identify certain generic 
features of the phonon spectra in the investigated fam- 
ily of Fe/Pt multilayers. To understand them better, we 
first summarize shortly the results of the electronic struc- 
ture calculations — they yield a tetragonal structure with 
a = b = 4.0246 A and c = 13.7906 A lattice constants. 
The former value is only slightly larger than iron lattice 
parameters (we found a — 4.0188 A for a test calculation 
with the same setup). This means that there are only small 
stresses in the multilayer — judging from the compress- 
ibility of the iron structure not higher than of the order 
of 0.3 GPa. The interlayer spacing results are the same as 
described above in section [21 

The equilibrium configuration was determined by com- 
paring the total energies of several possible stackings of Pt 
layer on Fe layer. The structure of the layers themselves 
was chosen following the experimental evidence that the 
bcc iron and fee Pt are present in this material. The crys- 
tal structure of the 4/4 system is given by the PA/nmm 
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Fig. T. Total phonon density of states p{uj) (solid lines) and 
partial phonon densities of states for Fe (dashed lines, red) 
and Ft (dashed-dotted, blue) atoms, respectively, as obtained 
for the FM ground state with: (a) 4/4 multilayer, and (b) 5/5 
multilayer, represented schematically in Fig. [T] 



space group, labelled as 129 (group Dj^ in Schoenflies 
notation), and for the 5/5 one it is PAmm, 99 (C|„) re- 
spectively. We have verified that the structure with the 
lowest energy is the one with a straightforward continua- 
tion of lattice points from bcc lattice of iron to fee lattice 
of platinum, as shown in Fig.[T]and described in section[21 

The calculated phonon dispersion relations for the 4/4 
and 5/5 structures are shown in Fig. [51 It is clear that the 
shape of the dispersion curves is quite similar, and the 
energy range of phonon spectra is the same for both sys- 
tems. The main difference is in the number of branches, 
due to the different symmetries of the investigated multi- 
layers. Thus one finds 3 acoustic and 27 (21) optic phonon 
branches for the 5/5 (4/4) multilayer, and some differences 
in the detailed dispersions of individual phonon modes. 
However, in Fig. [7] we can clearly see that these differ- 
ences do not affect significantly the distribution of spectral 
weight in the total phonon density of states p{oj), and the 
overall frequency ranges for the individual contributions to 
p{uj). This observation gives further support to the present 
choice of the supercell — even the a/2 x y/2 x 1 supercell 
used for the 5/5 system is large enough to obtain mean- 
ingful and representative results for the phonon density of 
states, while for detailed analysis of the dispersion curves 
it is necessary to consider instead a larger 2\/2 x 2 a/2 x 1 
supercell, as in the 4/4 system. 



The phonon density of states function p{uj) in Fig. [7] 
consists of two distinct contributions — heavier platinum 
atoms contribute with lower frequencies, while the vibra- 
tions of iron atoms are seen in the upper part of p{uj), well 
separated from the lower part. The structure of the upper 
band of the phonon density of states p{ijj) is clearly influ- 
enced by the contributions of the bulk-like vibrations of 
the inner iron layers, while the wide background is added 
by the interface atomic layers. 

While an experimental verification of the theory con- 
cerning the phonon dispersions in the Fe/Pt multilayers 
(Fig. El is not possible at the moment, we suggest that 
the projected phonon densities of states at iron atoms 
/9fc(<^) could be measured by the nuclear inelastic scat- 
tering, which make use of the Mossbauer effect [H]. Since 
the Mossbauer experimental techniques are sensitive only 
to iron atoms and depend on polarization for the consid- 
ered type of layered structure, we extracted the projected 
phonon densities of states for two different iron MLs in 
the 4/4 structure, with X — Y planar vibrations, and with 
perpendicular Z vibrations, see Fig. [8l As the projected 
phonon density of states for non-equivalent inner iron MLs 
of the 5/5 system are rather similar to each other, we 
analyse in detail /f)Fe(w) obtained for the 4/4 system, as 
it shows already the general trends representative of the 
considered structures. 

The differences between the X — Y and Z phonon 
polarizations in Fig. [5] are remarkable. While the spec- 
tral weight for the Z polarization is more uniformly dis- 
tributed, one finds characteristic maxima in the X — Y 
density of states near uj ~ 6.5 THz and uj ~ 8.8 THz for 
the central layer (bottom of Fig. [S]) . At the interface (top 
of Fig. [5]) an additional structure at w ~ 7.2 THz was 
identified, while the maximum at high energy is moved 
down to w ~ 8.0 THz, with low spectral weight above 8.3 
THz. Also for the Z polarization the spectral intensity in 
the high energy part oj > 8 THz is reduced. 

Altogether, the phonon spectra show the quasi-two- 
dimensional character of the investigated multilayer struc- 
ture, with large spectral weight in the range of 5.5 < w < 
8.5 THz for the interface, and a broader and extending up 
to ^ 9.3 THz projected phonon density of states pfc{^) 
for the inner Fe ML (Fig. |8]). This feature follows from 
the coupling to heavier Ft atoms at the interface, so the 
phonon frequencies are there reduced. The opposite situa- 
tion was found for Fe/FeSi interface [T4], where the spec- 
tral weight was shifted to higher frequency range for the 
modes which involve light Si atoms at the interface. How- 
ever, the projected phonon density of states found for the 
inner Fe ML in Fe/Pt system has several features (over- 
all distribution of spectral weight, the positions of max- 
ima) which make it similar to that of the inner Fe ML in 
Fe/FeSi multilayer. This confirms that the lattice vibra- 
tions depend predominantly on the local coordination of 
particular Fe atoms, which follows from the short-range 
character of interatomic forces. 
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Fig. 8. Partial phonon densities of states pFe(i^) projected at 
Fe atoms of the 4/4 Fe/Pt multilayer structure in FM config- 
uration for: interface (top), and central (bottom) Fe ML. In 
each panel the densities of states for X-Y and Z polarizations 
are plotted separately by solid (red) and dashed (blue) lines, 
respectively. 



5 Summary and Conclusions 

We have calculated the magnetic properties and lattice dy- 
namics of the set of Fe/Pt multilayer systems of varying 
thickness. In order to demonstrate the universal stability 
of the FM configuration between neighbouring Fe layers, 
we investigated spacer thickness dependence of the en- 
ergy difference between the AF and FM structures. This 
energy difference was quantified using the bilinear (Ji) 
and biquadratic (J2) coupling coefficients. It is quite re- 
markable that the values of Ji and J2 saturate fast with 
increasing thickness of Pt layer, indicating that the FM 
order of consecutive Fe layers is hardly influenced by the 
Pt spacer. These data may be therefore considered to be 
representative of larger systems as well. 

A systematic analysis of the magnetic structures shows 
that the magnetic properties of Fe/Pt multilayer are quite 
different from those of Fe/FeSi multilayers, where a crossover 
from the AF to FM interlayer coupling (between Fe layers) 
was found instead, with increasing thickness of the FeSi 
spacer layer p2]- In the case of Pt spacer the tendency 
towards FM order is well pronounced, and the magnetic 
moments of Fe atoms at the interface are enhanced. This 
agrees with earlier experimental findings which revealed 
rather robust tendency towards FM order in the Fe/Pt 
multilayers [5]. 



The phonon study has demonstrated that the phonon 
spectra consist of several modes, reflecting the geometry 
of a considered system with large unit cell. Platinum and 
iron vibrations are rather well separated from each other, 
and the vibrations of iron atoms contribute predominantly 
at higher energies to > 3.5 THz (w > 5 THz in case of the 
interface layer) , as shown by the respective partial phonon 
densities of states. The phonon spectra projected on iron 
atoms depend rather strongly on the multilayer geometry, 
and differ distinctly from those of the bulk Fe [19| , both 
for the interface and for the inner Fe MLs, at least when 
the iron layer consists of not more than flve atomic MLs, 
as considered in the present study. We believe that the 
maxima in the partial phonon density of states projected 
on iron atoms shown in Fig. |S] should be fairly easy to ob- 
serve and verify experimentally. Therefore, we hope that 
the present theoretical predictions concerning lattice dy- 
namics will help to understand the intricate dependence 
of the phonon spectra on the multilayer structure. 

Summarizing, the Fe/Pt multilayers exhibit not only 
interesting magnetic properties, but have also non-trivial 
and rich lattice dynamics, characterized by rather strong 
dependence of the partial phonon intensities on the po- 
sition of an Fe monolayer in the structure. It is believed 
that the presented results are representative for the Fe/Pt 
systems and will stimulate future experiments. 
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